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Abstract The generating function which counts partitions with the Plancherel 
measure (and its q-deformed version), can be rewritten as a matrix integral, which 
allows to compute its asymptotic expansion to all orders. There are applications in 
statistical physics of growing/melting crystals, T.A.S.E.P., and also in algebraic 
geometry. In particular we compute the Gromov-Witten invariants of the 
X p = 0(p — 2) © 0(—p) — > P 1 Calabi-Yau 3-fold, and we prove a conjecture of M. 
Marino, that the generating functions F g of Gromov-Witten invariants of X p , come 
from a matrix model, and are the symplectic invariants of the mirror spectral curve. 



1 Introduction: partitions 

Partitions are extremely useful in many areas of physics and mathematics. A partition 
with at most N rows, is an ordered sequence of N non-negative integers: 

A = (Ai > A 2 > ... > \ N > 0) (1.1) 

|A| = ^ Aj is called the weight of the partition A. The length of the partition n(A) is 
the number of strictly positive rows n(A) = #{Aj > 0}. 

It is also convenient to represent the partition rotated by 7r/4, i.e. we denote: 

hi = \i-i + N , hi > h 2 > . . . > h N > (1.2) 
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The hi's are in 1-1 correspondence with the downright edges of the rotated Ferrer 
diagram. 

Then, for various applications, ranging from statistical physics of growing/melting 
crystals [5T], T.A.S.E.P. [3SJ[T5], to algebraic geometry [45], one would like to measure 
partitions with the Plancherel measure V(X) (described below), and compute partition 
functions of the form: 

-k 

Z N (q;t k )= V W <?' A| e-v^^c.M (L3 ) 

n(X)<N 

where Cfc(A) are the Casimirs of the partition A (see section [L2|) . 

Here in this article, we find all the coefficients in the large q expansion of Z: 

oo 

In (Z^(q; tk)) = / J q l 9 F g (tk) + exponentially small (1.4) 

9=0 

Fq and Fx have been known for some time [17], and here, we prove that all the other 
orders Fg's with g > 2 are the symplectic invariants defined in |20j . 

Beyond computing the partition function, one may also be interested in finding 
the large size asymptotic shape of those random partitions. The shape of the typical 
partition has been known for some time [561 ESI EH EH E3 E3 S3 HZ] , it is related 
to Tracy Widom law [2S]- Here in this article, we recover the shape of the typical 
partition, and we compute corrections to all orders. We also compute all orders 
expansions of density correlation functions. 

We are also interested in a q— deformed version of the Plancherel measure V q (X), 
which has applications to algebraic geometry/topological string theory, in particular in 
computing the Gromov-Witten invariants of some Calabi-Yau 3-folds X p described in 
section 14.4.21 below (see [T2"| I4"T]). Here in this article, we prove, as conjectured by M. 
Marino |41] . that the Gromov-Witten invariants of X p are given by a matrix model, 
and more precisely by the symplectic invariants of [20] computed for a spectral curve 
of the form H(e x ,e y ) = 0. This is typically the form of a mirror spectral curve in the 
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type B topological string theory. Moreover, using the recent result of [IT], this also 
proves that this model is equivalent to the Kodaira Spencer field theory. 



Outline of the article: 

• Section 1 is an introduction, where we recall the definition of the Plancherel 
measure, the Casimirs. 

• In Section 2, we rewrite the partition function Z of the Plancherel measure, as 
a normal matrix integral. As an immediate consequence, we find that the topological 
expansion of In Z is given by the symplectic invariants of [20J , associated to a spectral 
curve which we compute explicitely. This also gives all orders corrections to density 
correlation functions. 

• In Section 3, we repeat the same derivation as in section 2, but for the g-deformed 
Plancherel measure partition function. We also rewrite it as a matrix integral, and as 
a consequence we find that its topological expansion, is again given by the symplectic 
invariants of [20J. In particular, we study the consequences for the Gromov-Witten 
invariants of the Calabi-Yau toric 3-folds X p . We prove the conjecture of [H] that the 
Gromov-Witten invariants, are the symplectic invariants of a mirror spectral curve. 

• In section 4, we explain the main applications: crystal growth, T.A.S.E.P., 
length of the increasing subsequences of a random permutation, and algebraic geometry 
Gromov-Witten invariants. 

• Section 5 is the conclusion. 

1.1 Plancherel measure 

The Plancherel measure is: 



where dim A is the dimension of the representation of the symmetric group £(|A|), 
indexed by the partition A. The Plancherel measure depends only on the partition A 
and it does not depend on N > n(A). 
Its q— deformed version is: 






where [h] is the g-number: 



[h] 



h/2 _ h/2 




i.e. 
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The Plancherel measure appears in many physical and mathematical problems. It 
is the natural measure on partitions. 

Among famous physical problems related to Plancherel measure are the 2D grow- 
ing/melting crystal [13133], the length of the longest increasing subsequence [S3], or the 
Totally asymmetric exclusion process (T.A.S.E.P.) [55], [To]. Plancherel measure also 
plays an important role in string theory and algebraic geometry, in the computation 
of Gromov- Witten invariants g H |4"6 | HTl W\\ . 

We develop those examples in section HI 
1.2 Expectation values and Casimirs 

In most applications, one would like to compute expectation values of the following 
form (as well as the g-deformed version of this): 

z N (g k )= e v ^ n e " Efc9fcCfc(A) ( lq ) 

n(A)<7V i 

where Cfc(A) is the k th Casimir, which can be defined as the k th term in the small z 
expansion of: 

1 N -(N--)z i 

E^^(A) = E<=*-" + '' + Wr-i (1 ' 10) 

k ' i=l 

The Casimirs depend only on A, they don't depend on N. 
We may write: 

- 1 

C k (\) = "£(hi ~(N- -)) k + C k>0 (N) (1.11) 



i=i 

For instance: 



n ; ^ 2 24 1 1 24 v ; 

i 

C 2 (A) = h i ~ ( 2N ~ !) E + \ N ( N ~ ~ !) = E X ^ + (1-13) 

i i i 

In general J2k 9kCk(X) can be written: 

Y,9kC k {X) = Y,~^Y, h " = E^) ( L14 ) 

k kit 

and it is important to notice that the coefficients t k may depend on N. 
Finally, what one would like to compute is a sum of the form: 

N d+1 

Z N oc V W U e ~ A{hl) ' A(h) = J2tkh k (1.15) 

n(X)<N i=l k=0 
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2 Plancherel measure 



In this section, we rewrite the sum 



1 — k 

Z N (q,h)= V(X) e-^iV-^c kW (216) 

n(X)<N 

as a matrix integral, and as a consequence, we obtain its topological large q expansion: 

oo 

In Z^tk)-^^ 1 F M ( 2 - 17 ) 

n=0 

2.1 Transformation into a matrix integral 

Consider the contour C which goes above the positive Real axis in the negative direction 
and below the positive real axis in the positive direction, i.e. it encloses all non negative 
integers. 



fa 



H 1 1 hH 1 1 1 1 1 h 

2 3 4 5... 



The function: 



ire-'** 



/(0 = -m-Om e-** = — Vtt (2-18) 

sm «) 



has simple poles with residue 1 at all integers (£ G Z). 
Therefore we have: 



-MM 



7<j i 



e ri( ft .-^) 2 n 



-A(hi) 



h\ 2 

hi,...,h N =0 i<j i 1 

m E rw^ n 



p -a(/k) 

/ii>...>hjv>0 



= iV! ^ P(A) ffe- A(hi) (2.19) 

n(\)<N i 

In other words, the Plancherel measure partition function is a normal matrix in- 
tegral: 



fei>...>hjv>0 ->H N (L) 



2 The ensemble Hjy(C) is the set of normal matrices whose eigenvalues belong to C, for instance 
iijv(M) = H n is the set of hermitean matrices. 
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with the potential: 

y/qV{x) = lnT( v /gx) - laT(-y/qx) + iir^/qx + In (y/qx) + A(y/qx) (2.21) 




And we write: 



A{y/qx) = ^fq - x In q + 



fc=2 



/.' V'' y/g 



^ 2 ^ +A 



y/qU{x) 



(2.22) 



Therefore, the derivative of the potential can be written (using Stirling's formula): 



V\x) 



ip(<Jqx) + ip(—<Jqx) + in H — + U'(x) 

y/qx 



2\n(^x) + U'(x) + — 



2n 



2 In x + 22 h+l 



(2.23) 



k=l 



y/qx ^ n(y/qx) 



2n 



X 



B 



2n 



y/qx ^ n(y/qx) 



2?i 



where S n are the Bernouilli numbers, and ?/> = T'/T. 



2.2 The rescaling factor and topological expansion 

We have introduced a rescaling factor ^fq such that hi = y/qxi, because typical h^s 
are not of order 1, they are large, and with this rescaling, we expect the typical x^'s to 
be of order 1 . 

Another reason for introducing a rescaling factor ^/q, is that it will serve as an 
expansion parameter. 

If the matrix model integral has a so-called "topological" large q expansion [541 [TO], 
i.e. an expansion in powers of q^ 1 here (and this must be the case by definition for 
formal series), then the coefficients are given by the solution of loop equations first 
found in [T8J, [13] , and are given by the symplectic invariants [20] . It was found in [20] , 
that the whole 1/q expansion can be written in terms of the spectral curve, i.e. the 
" large q density" : 



oo 



hiZ N ~2^Q 1 ~ aF B 

9=0 



(2.24) 
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where F g are the symplectic invariants of the spectral curve. Here, we shall see that 
the _F 5 's don't depend on N and on q. 

Therefore in what follows we aim at determining the spectral curve corresponding 
to the potential V(x) and to the integration contour C. 

2.3 The spectral curve 

The spectral curve is more or less the large q density of the h^s. It is related to the 
"large q limit" of the resolvent W{x) (see [20J ) : 



W(x) = "lim" — < Tr — - — >= "lim" / V 

yfq x — M \i 



, , v2.25) 
— hi 



Fixed filling fractions spectral curves can be found by the following recipeo 

we look for a set of cuts [a$, 6j] in the complex plane, and we look for a function 

W(x) analytical outside the cuts, which behaves like -^=^ at oo, and such that on each 

cut: 

W(x + iO) + W(x - iO) = V'(x) , Vxelauh] (2.26) 

and 

- — / W(x)dx = €t (2.27) 
2m JlaiM 

where e, are the filling fractions. Once the filling fractions are given, the solution is 
more or less unique. Therefore, we have to specify filling fractions. 

In fact, there are as many solutions of loop equations of matrix models (i.e. as many 
possible spectral curves corresponding to this potential), as the number of integration 
contours where the integral converges. A choice of integration contour, is equivalent 
to a choice of filling fractions. 

Here we want the integration contour to be C N , and we have to find which filling 
fractions it corresponds to. 

The spectral curve is the "large q" resolvent, i.e. it gives the asymptotic shape of 

a typical large partition. For this problem, the typical shape of large partitions has 

been determined for a long time [55l H6l Wf\ W2\ , and has a very special property called 

"arctic circle". We describe it below. 

3 The true definition of the spectral curve, is obtained by dropping the connected 2-point function 
from the loop equation, see [20 . But for 1-matrix models, this is more or less equivalent to solving 
the saddle-point equations for the his, and this simplifies to give the recipe presented here. 
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2.3.1 Arctic circle 

Our partition function Zjy = XwAKiV^'M e~ A ^ depends on N through the summa- 
tion bound. What was observed [52], is that typical partitions have a typical length 
< n(X) >= n(q, tk), and if we choose N > n, the h^s with N > i > n (i.e. hi < N — n) 
of a typical large partition are "frozen", i.e. hi = N — i, i.e. Aj = 0. 

n is sometimes called the arctic circle, because everything above the arctic circle 
is frozen [32J. This means that asymptotically Z N ~ Z w , up to exponentially small 
corrections, i.e. the F 5 's should not depend on N, and the N dependence should be 
only in the exponentially small corrections: 

oo 

In Z N ~ ^ 9F 9 + 0(e~^ c ^) (2.28) 

9=0 

This also means that, order by order in a 1/q expansion: 

In ^ ~lnZ w ~lnZ ?r ~ ^g 1-9 ^ (2.29) 

9 

In other words, if the arctic circle property holds, we may determine the F g J s by 
choosing N = V for instance, or by choosing our favorite value for N. We can also try 
to find a spectral curve for any N, and then check that indeed the _F g 's are independent 
of N, and this is what we do below. 

2.3.2 Chiral ansatz 

We try to guess the form of the spectral curve, and then we have to check the consis- 
tency of our hypothesis. 

The fact that all the hiS beyond n are frozen means that the large q resolvent (i.e. 
the spectral curve) can be written: 

TV 

W{x) = ± £ +W(x-b) , (2.30) 

V q i=n+l X "Vf" V? 

where W(x) has a cut [0, a] with an edge at 0. 

In general, W may have several cuts, depending on the tk's. Here we assume that 
the ife's are such that W has only one cut. We discuss this asumption below in section 

Thus, W(x) corresponds to a multicut solution, with N — n + 1 cuts, one of them 
is [b, b + a] with filling fraction -^=, and the others are degenerate cuts [^p, with 
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filling fraction -^=. This form of the spectral curve, is also called chiral ansatz in the 



litterature, after the solution of [THE?], see [12] . 

W(x) is thus a one cut solution and behaves for large x as 



~ n 3 

W(x) (2.31) 

and is such that 

2 - 1 

W(x + zO) + W(x - iO) = V'(x + b) V ^— , Wx e [0, a] (2.32) 

In other words, we have to solve the same problem as for W, but with a new condition 
that the cut is of the form [0, a]. 
Notice that: 

2 y> 1 - 2 \n X + b 1 | 1 f B2m ( 1 



y/q x + b — ^jjjr x yfqx y/q(x + b) m q m \x 2m (x + b) 2m 

(2.33) 

where B m are the Bernouilli numbers, and therefore the RHS of Eq. (I2.32p is: 



V'(x + b)- — 



2 A 1 



V i l=n +l y/q 

= 21111 + X>« I 1 - ^) + ^ - < 2 - 34 > 

Notice that it is independent of N. 
2.3.3 edge at zero 

We want to determine W(x), which has only one cut [0, a]. 

One cut solutions are better written in the Zhukovski parametrization: 

x(z) = j(z + 2 + -) , a = 4 7 (2.35) 

z 

Notice that each value of x corresponds to two values of z, namely z and 1/z. We call 
the domain \z\ > 1 the physical sheet, and \z\ < 1 the second sheet. 
Then we write: 

d / l v k d 

t k+ x ( x - n — i ) = + z- k ) (2.36) 

fc=l V / k=0 

where the coefficients u k s are functions of 7. We determine 7 by the condition: 

7 = e~ u ° (2.37) 
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The 1-cut resolvent W(x) is then: 

d CO j-y 

wW = *(,(,)) = + 21n (1 + 1 A ) + ^-^ - g 2n( ^. )) >, (2-38) 

indeed it satisfies Eq. (12.261) : 

d / — i \ k 1 oo „ 



and it has (order by order in 1/y/q) no singularity in the physical sheet \z\ > 1, and it 
behaves like 0(1/ z) for large z: 

y/qx(z) V. t(z) 

i.e. we find n: 

n-i = (wi + 2) 7 V? (2-41) 
The t^'s are then found from Eq. (12.360 : 

£ u fc («* + z- k ) = tk+i e~ kuo (z + -- Ul ) k (2.42) 

k=0 k=l 

In particular, the coefficient of z° and z 1 are: 



2«o = ELiC-i)***-!*-*" Efi? jSfin 

«i = -Eti(-l) fc Wie- feU0 E 



(2.43) 



i=0 j!(j+l)!(fc-2j-l)! 



which means that (after eliminating u\) Uq satisfies a transcendental equation of the 
form 

Polynomial (M ,e" M0 ) = (2.44) 

It is important to notice that this equation is independent of iV and q, and thus, uq, U\, 
and all the u^s depend only on the t^s, they don't depend on N and q. 

The resolvent W(x) is thus given by Eq. (12.301) : 

l JL i l 
W(x(z)) = — — + , =b + 1 (z + 2 + -) (2.45) 
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The spectral curve is the density y = ^(W(x— iO)— W(x+iO)), i.e. the discontinuity 
of W along the cut, i.e. y(z) = \{ui{l/ ' z) — uj(z)), therefore our spectral curve is: 




(2.46) 

The same spectral curve was obtained by Marshakov and Nekrasov in [32]) it was 
derived from finding the optimal profile function which maximizes the free energy Fq. 
It turns out that the extremization equation gives also Eq. (12. 26ft . 



2.3.4 Consistency of the chiral ansatz 

We have to check that the F 5 's are independent of N. The resolvent W(x) is paramet- 
rically given by the spectral curve: 



x(z) 



N 



[z + --u x ) 

z 



+ 1 y/qx(z)-N+i 

Eoo B 2n 
71=1 



+ ELi^- fc + 21n(l + l/^) + 



(2.47) 



2^qx(z) 



It was proved in [20], [21], that the _F 9 's are invariant under transformations of the 
spectral curve which preserve the symplectic form dx A dy (this is why they are called 
symplectic invariants). In particular they are invariant in we add a constant to the 
function x(z), and if we add to W(x(z)), any rational function of x(z). Therefore, the 
previous spectral curve is symplectically equivalent to: 



x[z\ 



e- U0 (z + 



'2.48) 



and this last spectral curve is clearly independent of N and q. 

Therefore the _F 9 's are indeed independent of N, which proves that the chiral ansatz 
was consistent. 

As a bonus, we also find that the F 5 's are independent of q. 
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2.4 Summary: asymptotic expansion 

Using the known properties of matrix models, the logarithm of the following partition 
function 

i-fc 

Z N (q,t k )= V(X) q W e~^ tj ^ c ^ (2.49) 

n(X)<N 

can be written as a large q expansion: 



lnZ N (q,t k ) ^J^Q 1 ^ F n(£{tk)) 



n=0 



(2.50) 

where the _F n 's are the symplectic invariants of [20] computed for the spectral curve 
£(tk) of Eq. f)2.46p . The specral curve, which depends only on the ti-'s, is found by 
solving the transcendental equations Eq. ( 12.431) . 

Some explicit examples are treated below in section 12.61 



2.5 Asymptotic shape of large partitions 
2.5.1 Leading order shape of large partitions 

If we write z = e^ we have the equilibrium density in trigonometric form: 



" lim " ^E(^-^ = - = -(^+E^ sin (^) 



x = & + 2 7 (l + cos0) (2.51) 
and the integrated density I = ^fq J^ +Al p e qdx is: 



^2(sin0 — 0cos0) + Ui{4> — sin cos < 



7T 

d 



k=2 



k + 1 k — 1 



— I 

2" 



so that 7(0 = 0) = and 1(0 = tt) = n 

I(x) is the integral of the density, i.e. the number of h^s between x and b + 4 7 . 
Remember that hi = ^/qXi, and thus, hj^ = ^fqx = ^/q(b + 2 7 (1 + cos0)). In 
other words, the inverse function ^fqx(I) gives the typical partition hi( x ). We have 
parametrically 

hi — N — n + 2 7v /g(l + cos 0) 



' COS 0) + Ml (0 - Sin COS 0) - J2k=2 U k (j 



(2.53) 



/ _ [ 2 f Sill v.-/ v. v ■ v.-/ , " ! ' ' > V/./ V w, ■ v.,/ / / .->"/,' i /._,_] / . _ | 
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The typical Xi = hj + I — iV's are given by: 



\ I = 7-71 + 27^(1 + 008 0) (2.54) 

If we plot Xi + I as a function of A/ — /, we get the typical shape of the 7r/4 rotated 
partition. Examples are plotted below in section 12.61 



2.5.2 Subleading fluctuations 

The correlation functions wlf^s defined in [20] give all the corrections to the densities 
and density correlations of the h^s, to all orders in ^fq. For instance we have: 

oo 

<E^37r>-E^^ (9) w ( 2 - 55 ) 

and more generally the cumulant functions expand as: 

nf^T^)) -Efl^W?^,...,^) (2.56) 

V fc=l \ i Xk VV / / c 9=0 

where wfi s are the correlators defined in [20] for the spectral curve Eq. ( 12. 46ft . 

By taking the discontinuity (or by substracting the principal part), we also obtain 
the correlations of densities p(x) = J2i $( x ~ ^«) : 



( \ i 

ITT /.ITT < ' 

and more generally 



in 2in 

9=1 



(p( Xl ) . . . p(x n )) c ~ E 9 1 "*" 5 ^(^i. • • • , *n) (2-58) 

(2V2 7rr 9= i 



The method for computing the Wn s is explained in [20], and it can easily be imple- 
mented on a computer. 



2.6 Examples 
2.6.1 No Casimirs 

We choose all = 0. 

Consider the partition generating function: 

Z N (q)= E nA)? |A| = l+? + y + ^ + --- (2-59) 

n(X)<N 
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Equations Eq. (12.431) reduce to Uk = for all k, and 7 = 1, and n — | = lyfq. The 
spectral curve is thus: 



and therefore: 



x(z) = 2 — yH +Z + - + 2 

£= J 

y(z) = \n(z) 



]nZ N (q)~J2 <l l ~ 9 F 9 {S) 

9=0 



(2.60) 



(2.61) 



In fact, it is well known [46J that when iV — >• oo, the generating function of partitions 
is: 

^oo = ^P(A)g |A| =e" (2.62) 



i.e. 



which shows that: 



In Zoo = q 
F 9 (£) = $ 9 ,o 



(2.63) 
(2.64) 



Density and typical shape 

We have: 

/ = — — (sin — cos ( 

7T 



A = I - - + 2 v /gcos( 



(2.65) 



i.e. by writing A — I = l^fqu we recover the famous [HI [56] asymptotic shape of large 
partitions: 



X-I = 2 v ^u , A + / = Ay/q ( v 1 — w 2 + wArcsin('u)) (2.66) 



A+I 




A-I 
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2.6.2 2nd Casimir 



In this example, we choose only t 2 7^ 0. Consider the following sum: 



n(X)<N 



Equations Eq. (I2.43P give: 

2u = -t 2 e~ u ° m 
i.e. mq is determined by: 



Mi = t 2 e 



-Uq 



[2.67) 



(2.68) 



- 2m e 2u ° = t\ 



- 2u ° = J2^r t i k = t i + 4 + ^4 + --- (2.69) 



fc=l 



Then we have: 

Mi = V-2m , 7 = e" 
The spectral curve is: 



n ~2 



(2 + Mi)e-"° y/q (2.70) 



S(t 2 



x(z) = e~ u ° (z + z' 1 - v^mo") 
V^2^ ( 1 



y(z) = \nz+ - 12 



2m = t\ e~ 2u ° 



(2.71) 



The same spectral curve was already found in [22] who computed F . Here, we also 
have the full topological expansion: 



\nZ N {q,t 2 ) r^Y,^' 9 F g (S(t 2 )) 

9=0 

For example we find (using mathematica with the definitions of [20J): 

Fi = ^ln(e- 2 "°(l + 2M )) 
e 2u ° ul (1 - 12m ) 
2 180 (1 + 2m ) 5 

As an example, we plot the typical shape of a partition with t 2 = 0.2: 

A + I 



(2.72) 



(2.73) 




x-i 
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and corrections to the limit shape are given by: 




y/qx(z) - hi 



u(z) + -^\z)+0{q- 2 ) 



(2.74) 



where 




2.7 Phase transitions and number of cuts 



In section 12.3.21 we made the assumption that W(x) has only one cut [0,a\. Then we 



could compute the density of the h^s on the cut. It may happen that the density we 
compute this way is not positive, which means that our 1-cut assumption was possibly 
wrong. 

For example, in the example of section 12.6.21 a "pure gravity" phase transition (i.e. 
conformal field theory (3,2), characterized by the Painleve I equation) occurs when 
t 2 = e" 1 (i.e. u = -|). 

In fact, all situations may occur as we vary the parameters t^'s: cuts may split, 
merge, shrink, or appear, they can move to the complex plane, and branch so that 
they form trees. Most of those transitions have been described in the random matrix 
theory literature, and take place also in this particular matrix model. 

In principle, the correct assumption, is the one which minimizes ReFo. 

It would be interesting to study the phase transitions between different regimes of 
the tfc's, and see how the cuts change, and thus how the shape of typical large partitions 
change. This would be merely an adaptation of mostly known phase transitions of 
matrix models, specialized to this case. We won't do it in this article. 

However, a non-positive density is not a problem when the partition functions we 
consider are only formal series. 

The importance of the 1-cut assumption, comes from the fact, that most often, 
this is the assumption which gives the formal generating series which enumerates some 
objects. Formal generating series are very often obtained from Feynman's method, 
i.e. a perturbative expansion near a minimum of an action. We keep the quadratic 
part of the action in the exponential, and expand the non-quadratic terms, so that we 
end up having to compute series of polynomial moments of Gaussian integrals, which 
are represented diagrammatically through Wick's theorem. This method naturally 
associates combinatorics objects to integrals [19]. 
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A 1-cut case means that all eigenvalues of the random matrix are expanded near 
the same minimum, whereas multicut cases mean that we expand near several minima. 
Multicut cases lead to formal series which enumerate objects of different types (one 
type for each minimum, i.e. for each cut). This is less natural than 1-cut cases where 
we enumerate only objects of the same type. 

Thus, for formal series, the correct assumption on the number of cuts, is not ob- 
tained by minimizing ReF , but is obtained from the type of objects which we wish to 
enumerate, i.e. near which minimum we want to expand. 

Here in this article, we focused on the 1-cut case for simplicity, but it is clear that 
the entire method would work for any number of cuts. The symplectic invariants of 
[20] were indeed defined for any number of cuts. 

2.7.1 Universal regimes 

Since we have a matrix model, we find the same universal regimes as usual in matrix 
models. Correlation functions in various regimes are obtained from some universal 
kernels. Also, it was found in [20] (cf theorem 8.1 in [20]) how to study the vicinity of 
phase transitions of the F g 's. 

For instance it is clear that we have "sine-kernel" in the bulk of the spectrum at 
microscopic scaling, "Szego kernel" in the bulk at macroscopic scaling, "Airy kernel" 
and Tracy- Widom law [58] (which is also the conformal minimal theory (1,2)) near 
regular turning points, and we have the Kadamtsev-Petviashvili (KP) kernels (i.e. 
minimal conformal field theory (p,q)) near singular points of the type y ~ x p ^ q , and 
polynomial kernels of [23] when a new cut opens. In general, the KP kernel for a (p, q) 
critical point, is the kernel (defined in section 9.1. in [20J) of the limit spectral curve 
(see section 10.3 of [20]), of equation 



However, we shall not study the details of those universal regimes and transitions 
in this article. 



3 The q-deformed Plancherel measure 

In this section, we repeat the same transformation into a matrix integral, chiral ansatz 
and topological expansion, for the g-deformed Plancherel measure. 



'p,<i 



' = Q(z) 
y(z) = P{z) 



deg Q = q 
deg P = p 



(2.76) 
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We compute: 



k 

Z N (q,t k )= e-^^S^ Cfc(A) q = e~*- (3.77) 

n(X)<N 

in the small g s topological expansion regime: 

oo 

hyZ N (e- B ',t k )~STg%>-* F g (3.78) 

3=0 

We shall discuss an important application of this model to topological string theory 
in section 14.4.21 below. 

3.1 Transformation into a matrix integral 

Let 

q = e-» , T = g s (N-~) (3.79) 



Let the contour C q be a circle of radius r with \q l \ > r > 1 centered in 0. In 
particular it encloses all the q h with h G N. 




Consider the g-series: 

oo 1 

g(x) = fj(l " ~q n ) (3-80) 

n=l 

It vanishes when x = q h with h G N*, and we have: 



/i\2 iwh -h(h-l)/2 

and thus the function: 

«(l) 2 e-f lnW e^ (ln " )2 
(1 - x) ^/xg{x)g{l/x) 

has simple poles with residue 1 when x = q h for h G N (one easily sees that it has a 
pole of residue 1 at x = 1 as well). 
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Notice that if x = q h with h G N, we have: 

h 



[h]l = H(q-^ - q^) = q -«>+M (3.83) 
Therefore we have: 



-^(ln^^/lnq) ffl 1 /^) p 21^ 



^ r<7 t - r > 1 



-A(fti) 



hi,...,h]s/=0 i<j i L lJ 

hi>...>h N >0 l l >' i 



h i~ h j h j~ h i s n 

2 > z 



(3.84) 



e n - <j " - n^ (t " 

/ii>...>/ijv>0 L * J ' i 



In other words, the g-deformed Plancherel measure partition function is a matrix 
integral: 



T Ui<i lt„ kj]2 T]e- A(hi) = ^ f dMe-^ v(M) (3.85) 



with the potential: 

— V(x) = \ng(x) — IngCi/x) + (N — l)hix -\ lnx + ln(l — x) 

9s 9s 

+- Y (~ l f tk (T + lnx) fc (3.86) 

Remark: since the contour C q is a circle of radius r, The ensemble H N {C q ) is (up to 
a scaling r) the ensemble of unitary matrices: 

H N (Q=rU(N) (3.87) 

It was already noticed [5J [33] , that sums over partitions with the Plancherel measure 
have many common properties with unitary matrix models, and here the connection 
becomes very apparent. 
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3.2 The potential 

We use the following asymptotic expansion of the g-series (for > \x\ > 1): 



1 00 R 
In (g(x)) = ~-J2 L^-mix' 1 ) 9? (3.88) 



q s * — ' ml 

ys m=0 



where B rn is the m th Bernoulli number (B = 1,B\ = — |, . . ., and .£>2m+i = if 
m > 1), and Li n (x) = ^2^ =1 x k /k n ^ s the polylogarithm function. We have: 



g'(x) 
x- 



1 °° R 

- ^ Li^x- 1 ) ^ ^ (3-89) 



The polylogarithm functions have many properties, in particular: 

Li' n = -Li n _i (3.90) 

Li-m(l/x) = (-l) m Li_ m (x) m>0 (3.91) 

Li^x) = -ln(l -x) , Lio(ar) = — — (3.92) 

X Ob 

and Li_ m with m > is a rational fraction with poles at x — 1 and which vanishes at 
x = and at x = oo. 
Therefore: 

-ins B Tt 



xV'(x) = y^(Li 1 ^ m (x)+Li 1 ^ m (x 1 ))—^g 



m=0 

d 

+ (iV - l)g s + iir + g s -^j - ^(-l) fe t k+1 (T + lux)' 

k=0 

(-21n(l-x)+lnx-i7r) + 

oo n 

+2 5>i- am (*)|^ 

m=l 

d 

+ (iV - 1)<? S + Z7T + ^— - ^(-l) fe t k+1 (T + lux)' 

Ob X 

k=0 

d 

-2 In (1 - x) + In x + T - ^(-l) fe 4+i( T + lnx) fc 

k=0 

oo „ 



I + 2j]Li 1 _ 2m (x)^^ 



(3.93) 



Notice that terms on the last line are rational fractions with poles at x — 1, and which 
vanish at x = and at x = oo. 
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3.3 Restriction to 2 Casimirs 

In principle, the recipe above, should allow to find the spectral curve for any tk's, but 
for simplicity, from now on, we shall assume that tk — if k > 3, and we write: 

ti = t , h=p-l , t 3 =U = t 5 = ... = (3.94) 

The potential is then: 

xV\x) = -21n(l -x) + p{\nx + T) -t 



(3.95) 



2m! s x 

m=l 



Those notations (t,p — 1) = (tx,t 2 ) are motivated by the application to algebraic 
geometry [12j UT], namely the computation of the Gromov-Witten invariants of X p , 
which we present in section 14.4.21 see eq J4.135[ 

Also, a motivation for the assumption = for k > 3, is that the spectral curve 
is much simpler in that case. For higher t^s, we would need to play with tedious 
identities among polylogarithms, and thanks to this assumption, we don't need higher 
polylogarithms than Li\{x) = — In (1 — x). 

3.4 Chiral ansatz and dependence on N 

We shall not repeat the discussion of sections 12.3.11 12.3.21 and 12.3.41 We let it to the 
reader (following the same steps as in the undeformed Plancherel measure) to check 
that the spectral curve has an arctic circle property, i.e. all h^s with i > n are frozen, 
and we can choose a chiral ansatz [HI [12] such that the F 9 's don't depend on N. 
Therefore we chose N = n, such that the edge of the small q average distributions of 
the hiS, is at h = 0, i.e. x = 1. The authors of [12] checked the consistency of the 
chiral ansatz, by verifying that Fq is independent of T (there is a T 2 scaling between 
ours and that of [T2]). 

In other words, in an appropriate regime of t and p, we may choose T, such that 
the spectral curve corresponds to a 1-cut distribution, with edge at x = 1. 

3.5 Spectral curve 



Like in section I2.3[ the spectral curve is determined by: 

W[x + iO) + W(x - iO) = V\x) , V x e [a, 1] (3.96) 
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and W(x) is analytical outside the cut [a, 1], and at large x: 

W(x) ~ ^ (3.97) 

x 

Since we look for a 1-cut solution, we again use the Zhukovski parametrization and 
write: 

x(z) = 1-7(2 + 2 + -) , o = l-47 (3.98) 

Notice that each value of x corresponds to two values of z, namely z and 1/z. We call 
the domain \z\ > 1 the physical sheet, and \z\ < 1 the second sheet. 
Define zo such that \zo\ > 1 and x(zo) = 0, that is: 

I = (l + z )(l + l/z ) (3.99) 

7 

Then, the resolvent is W(x(z)) = u>(z): 

x(z)u(z) = -2(ln(l + l/z) -ln(l + l/z )) 

+pQn(l-l/^ )-ln(l-l/, 2o 2 )) 

+ ^U 1 - 2m (x(z)) — g s + - 

m=l \ \ / / 

(3.100) 

Indeed, one easily checks that ui(z) has no singularity in the physical sheet \z\ > 1. 
The condition u(z) ~ -%r at large z implies: 

2 In (1 + 1/zo) - pin (1 - l/z 2 ) = T = g s (n - ~) (3.101) 

Then, the loop equation Eq. (I3.96P implies: 

x(z)(u(z)+u(l/z)-V'(x(z))) = 

2(ln 7 + 21n(l + l/2 )) 
-p{\n( 1 z ) + 2\n(l-l/z 2 )) 
+t-pT 

(3.102) 

therefore z , T and 7 are determined by the 3 equations: 

= 2 (In 7 + 2 In (1 + l/z )) - p (In (jz ) + 2 In (1 - I/2 2 )) + t - pT 
T = 2 In (1 + l/*b) - pin (1 - I/2 2 ) (3.103) 

1 = (l + z )(l + l/^o) 
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and after a little bit of algebraic manipulations, zq is determined as a function of t and 
P b y : 




(3.104) 

One should notice that this equation is symmetric in p <-> 2 — p, which reflects the 
symmetry of the Calabi-Yau space X p «-> X 2 _ p . The same equation was found in [T2] 
where the authors wrote w = 1 — \ (see equation 4.34 in [12] or equation 3.14 of [4"T]). 
Order by order in e~* we have (using Lagrange inversion formula): 



k=2 



k\ 



n - 2 ) 

j'=o 



= e~* + p(p - 2)e" 2 ' + 

We also determine T and 7 through: 

(1 - -) p 



p(p-2)(3p(p-2) + l) 



-3t 



+ 



'1 + -i-)2-P 

20 / 



1 

7 



(l + z )(l + l/^ ) 



(3.105) 



(3.106) 



The arctic circle is at N = n 
Finally, our spectral curve is: 



2 T 9s- 



£(t,p) 



x(z) 



y{z) 



: z o' zz o' 
) 2 



— (1 - JlV(p-2) 
2 v _2' 



2, 











(3.107) 

One may check that this spectral curve is identical to that of [121 E] ( see ec L- 3.31 
in [JT]), therefore it gives the same F g 's. 
This proves the conjecture of [41J. 



3.6 Topological expansion 
3.6.1 Topological expansion 

Since Z is by definition a formal power series in g s (it is the generating function of 
Gromov-Witten invariants of X p , see section 14.4.21) . and is given by a matrix model, 
using [20], we find that Z has a topological expansion 

00 

\nZ = Y,9 2 s 9 ~ 2 F g {£) (3-108) 
9=0 
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where the coefficients F g {£) are the symplectic invariants defined in [20J. 

Let us compute the few first orders (the first orders were computed in |12j): 
• Fq is found from the general theory of [20] through its 3rd derivative: 

}2 Ei 1 _ 1 \2\ 



1, T{k{p-lY) th 



<9t 2 zg' ^ k\T(kp(p-2) + 1) 

It can be checked that this leads to the expression found in 



/•;, . /J ( / ,-2)/./: 3 ( r ^) + (p-i) 2 L»3(i/^)- p(p 2) i P 1)2 (log (1- 1)) 3 (3.U0) 
1 - 4 6 4 



For instance if p = or p = 2 we have Fq" = e kt , i.e. F = J2 k = Li 3 (e*). And 
for p = 1 we have Fq" = ^ k {—^) k Q kt , i.e. F = Lz 3 (— e*). 
• Fi is equal to (according to [20]): 

F t 4 in (TVdM-D) - ± in ( " (( ^-"' ) (3.111) 



which is the same as in [41] . 

• i<2 is obtained from [20] (using a few minutes basic Mathematica programm) 

F 2 = , ( (p ~ 1) 8 (-1 + 124 - 124) 

2880 {zl - (p - l) 2 ) 5 V F J K ° ; 

+(p - l)% 2 (-5 + 4 + 2^) + 35(p - 1)% 4 (-1 + 4) 

+(p - 1) 2 4(2 + 4 - 54) + 4(i2 - 124 + 4)) 

(3.112) 

• Higher F ff 's can be obtained easily with a higher computing time... 
3.6.2 Comparison with the conjecture of | |41| 



The authors of [12] have computed the spectral curve as the solution of the large TV 
saddle point equation for the h^s. Then, Marcos Marino [UJ computed the first few 
orders F Q and F\ for that spectral curve, and conjectured that the Gromov-Witten 
invariants were generated by the symplectic invariants F g of the same spectral curve 
for all g. However, our approach is rather different, and it is worth to comment on 
"why it works ?". It is due to a few miracles. 

In fact, in their computation, the authors of [12J, worked in the leading large N 
limit, and they made several approximations. 

- The first of them, is that they set T = Ng s instead of T = (N — \)g s . This makes 
no difference to leading order, but it can make a big difference to subleading orders. 



24 



- The second of them, is that they neglected all terms with a 1/N behaviour in the 
Plancherel measure, keeping only L12 and Li\, they neglected all the Bernouilli series. 

- The third of them, is that they neglected the discreteness of the h^s, and replaced 
them by continuous variables, i.e. they replaced sums by integrals, whereas here in 
this article, we encoded the discreteness of the hiS by the q— series g(q h ), which forces 
(without approximation) the h^s to be integers. A miracle is that the ^-function 
combines very well with the Plancherel measure. 

Out of those approximations, they deduced a spectral curve, and quite miraculously, 
it turns out that it is indeed the same curve as we found here, modulo symplectomor- 
phisms. In some sense the 3 approximations compensate each other, and they indeed 
found the right spectral curve which gives the answer to all orders. This miracle, is of 
course due to the very special nature of Plancherel measure, it would not have worked 
with another arbitrary measure. 

3.7 Symplectic invariance and mirror symmetry 

Here, we use the fact that the _F 9 's are symplectic invariants, i.e. two spectral curves 
which have the same symplectic form dx A dy have the same Fg's. 

Notice that the spectral curve Eq. (13.1071) has the form: 

f x(z) = u(z) 

£{t ' p) = { v(') = i W) (3 ' 113) 

where both u and v are rational fractions of z. 
The symplectic form dx A dy is thus: 

dx A dy = d\nu(z) A dkiv(z) (3.114) 

therefore the spectral curve £{t,p) is symplectomorphic to the following new spectral 
curve: 

= lnu(z) . u(z) =720(1- ^, 



8 <t>p) = S M 1 m n i/'-iV ( 3 - 115 ) 

y(z) = hiv(z) , v(z) = - 



which is an algebraic curve in e x and e y , i.e. there exists a polynomial H(u,v) such 
that: 

H{e x ,e y ) = (3.116) 
This is typically an equation of a mirror Calabi-Yau 3-fold in type B topological string 
theory [30l Efl EU [9] . 

Remark: notice that the transformation (u,v) — ► {w ,v) is also a symplectic trans- 
formation which leaves the _F ff 's unchanged. It corresponds to the so-called framing 
transformation. 
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3.8 Mirror curve 

The curve H(e x ,e y ) = can be written by eliminating z: 
(e* + e^)(l + Vv^ = W^) tl)! r ._ l(( i + Zo)( i + I )( i_e-)-2) (3.117) 

where Tj is the j th Tchebychev polynomial, such that Tj(z + z^ 1 ) = z j + z~ j . 

• Example p = (z = e^ 2 ): 

v + - + 2 = (l-u)(l + e*)(l + e = £) (3.118) 
v 

i.e. 

= (1 + e y )(l + e~ y ) + (e x - 1)(1 + e*)(l + e^) (3.119) 

• Example p = 1 (,2q — 1 = e*): 

+ ^ = (1 + ^o) 2 m + (1 - zofu- 1 - 2z 2 

• Example p = 2 (z = — e*/ 2 ): 

v + v- 1 + 2 = {1- u -1 )(l + e*)(l + e T) 

which is symplectically equivalent to the case p = by changing u 

• Example p = -1 {z 2 (l - z 2 ) 3 = e~*): 

v 2 + v - 2 = (1 + I)\ + (i _ l)2 u -i _ 1_ 

Zq Zq Z q 

4 Some applications 

In this section we discuss some applications of our method. 

Important applications of sums of partitions with the Plancherel measure include 
growth phenomenon (crystals), length of increasing subsequences of a random sequence, 
Totally asymmetric exclusion processes, and algebraic geometry through the computa- 
tion of Gromov-Witten invariants and Hurwitz numbers. 

Many of those applications are based on standard Young tableaux. 

A standard Young tableaux , is a Ferrer diagram filled with numbers such that 
the numbers are increasing along each row and along each column. The number of 
standard Young tableaux of shape A, is dim(A). 



(3.120) 



(3.121) 



u 



-l 



(3.122) 
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4.1 Longest increasing subsequence 

Let cr G Sfc be a permutation of k elements. Call 1(a) the length of the longest 
increasing subsequence of a. 

The Robinson-Schensted algorithm [S21 E3] gives a bijection between the group of 
permutations Sk and the set of pairs of Standard Young Tableaux with k boxes (T, T') 
with the same shape A (shape(T) = shape(T') = A, |A| = k). This bijection illustrates 
the formula: 

k\= (dim(A)) 2 (4.123) 

A,|A|=fc 

The length of the longest increasing subsequence is in bijection with the length of the 
partition n(X). The uniform probability law on permutations induces the Plancherel 
law on partitions. Therefore, the probability law of the length of the longest increasing 
subsequence of a random permutation (uniformly chosen), is the probability law of 
n{\) with the Plancherel measure [H El [33l [3H EH EH1 EH [53l |55l [56]. What one 
needs to compute is the expectation value: 

E(l(a) < N) = E fciiV = I J2 1 = k] E V ^ ( 4124 ) 

' creS fc ,l{cr)<N X, \X\=k,n(X)<N 

Instead of working with permutations of a given number of elements k = |A|, it may 
be more convenient to work with a " Poissonized" version (grand canonical ensemble), 
where we sum over all numbers k = |A|, with a "chemical potential" of the form g' A L 
In the end, we need to compute: 

k 

Z N (q) = e-«J2 E *,Ny_ = e ~* E r ^1 lM ( 4 - 125 ) 

k '' X,n(X)<N 

Baik, Deift and Johansson [6] proved that Z^iq) is related to a Haenkel determinant, 
and that in the large q limit, the distribution of the length of the longest increasing 
subsequence converges towards the Tracy- Widom law [58]. In particular they could 
find the "typical shape" of a large partition. For instance they found that the typical 
size k — |A| is given by: 

<|A|>=g (4.126) 
and the typical length 1(a) = n(A) is given by: 

< l( a ) >=< n(X) >= 2^q~ (4.127) 
and typical fluctuations are of order g 1 / 6 . 

< 1(a) - 1^ >~ 0(q 1/6 ) (4.128) 

Here, in this article, we compute not only the large q leading order of Z^(q), but the 
whole 1/q expansion of Z^(q) for large q. 
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4.2 Growing/melting crystal 

The Robinson-Schensted algorithm also gives a simple growth model for a 2- 
dimensional crystal. 

The crystal is represented by 7r/4 the rotated partition A. At time t G N, we drop 
a new box, falling from the sky at a random position a(t). When the box hits the 
crystal, it slides to the left or to the right until it can no longer move. If we assign the 
box, the time t at which it arrives (resp. the position a{t) from which it came), we 
clearly obtain at each unit of time a standard Young tableau, and the number of ways 
of obtaining a given standard tableau of shape A after | A| units of time, is the number of 
permutations a which can lead to it, i.e. it is the number of standard Young tableaux 
of shape A, i.e. it is dim(A). The 2 standard tableaux we obtain from assigning the 
box t or a{t) are precisely the 2 tableaux of the Robinson-Schensted bijection. 



This is a model of a growing crystal. 

The other way around, it is also a melting crystal, where at each unit of time only 
boxes which are at corners can evaporate. If we label each box by the time at which it 
evaporates, we get a standard tableau. 

Therefore it is natural to associate to each partition, a weight which is its entropy, 
namely the number of ways of obtaining it, i.e. dim(A) 2 . If the number of boxes is not 
fixed, we may consider a grand-cannonical ensemble, with a weight depending on A, 
for instance of the form g' A L Again, the partition function we need to compute is of 
the form: 



A 

and we are also interested in expectation values of various observables, typically ex- 
pectation values of Casimirs. 

This analogy between partitions and crystals has been very useful and fruitful, and 





(4.129) 
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has generated a considerable amount of works and discoveries in physics and mathe- 
matics [5H Eg El [251 EH HE]. 

4.3 T.A.S.E.P. 

The acronym T.A.S.E.P. stands for totally asymmetric exclusion process. It is a famous 
model of statistical physics, where particles are at integer positions on the real axis 
[381 [151 HS1 [29]. At each unit of time, we choose one particle, and with a certain 
probability, either it stays at the same place, or it jumps to the next position to the 
right if unoccupied. This model is very important because it is the simplest "out of 
equilibrium" statistical physics model. 

There is a bijective mapping between this model and the growing crystal. Indeed, 
in the growing crystal, at each time t, we have a partition A(i) of weight |A(t)| = t. 
If we consider the 7r/4 rotated partition, we have for each t, a set of ordered integers 
hi(t) > h,2(t) > ... > hpfit) > 0. If we interpret the integer hi(t) as the positions of 
the i th particle at time t, we see that those particles follow a T.A.S.E.P., whose initial 
condition at t = is a "Dirac see" (Z_ full of particles, N empty). 



h 5 h 2 h x 
t=0 ••••••• 




Here, in this article, we give a method to compute the correlations of particle 
positions of the form: 

((fi^)(fi^)"(fii^)> (4130) 

for all k, and to any order in the large q expansion. 

4.4 Algebraic geometry 
4.4.1 Gromov-Witten theory of P 1 

P 1 is the complex projective plane, i.e. the complex plane with a point at oo, it is (the 
only) compact Riemann surface of genus 0, it is also called the Riemann sphere. 
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A stable map (C, pi, . . . , p n , f) to P 1 , is the data of a (possibly stable nodal) curve C 
(of complex dimension 1), with n smooth marked points Pi, ■ ■ ■ ,p n , and a holomorphic 
function / of degree d, from C to P 1 . A nodal curve, is a possibly degenerate compact 
Riemann surface (some cycles may have been pinched), and stability means that if we 
remove the pinched cycles, all connected domains have strictly negative Euler charac- 
teristics, so that there are only finitely many holomorphic automorphisms preserving 
the marked points. The degree d of map / : C — > P 1 , is the number of preimages of a 
generic point in P 1 . Critical points are the points with less than d preimages, i.e. such 
that df = 0. Generically, if C is a smooth curve, / has 2d + 2g — 2 critical points. 

The set of all stable maps (C,p\, . . . ,p n , f) to P 1 , is a finite dimensional complex 
manifold (it is locally described by a finite number of complex parameters called moduli, 
indeed the data of a given degree function is more or less the data of the coefficients 
of a polynomial), which is called the moduli space. 

Let 7W 5in (P 1 ,c?) be the moduli space of stable maps of degree d to P 1 . Since each 
point Pi is smooth, we have a natural line bundle Ci over Aig^F 1 ^), whose fiber is, 
for each point (C,pi, . . . ,p n , f) G A / i 9i „(P 1 , d), the cotangent space of C at p i7 i.e. T*.C. 
We can consider the first Chern class c\(Ci) of this line bundle, i.e. the curvature form 
of a U(l) connection. 

There is a virtual fundamental homology class on A^ 9j „(P 1 , d), of dimension 2d + 
2g — 2 + n where the intersection theory can be computed (i.e. where classes c\(Ci) can 
be integrated). By Poincarre duality, classes are dual to cycles, and integrals of classes 
compute the number of intersection points of the cycles, provided that this number is 
finite, i.e. provided that the sum of dimensions is the total dimension of the whole 
space. If the dimension is wrong we say that the intersection number is zero. The 
virtual class, allows to extend the generic dimension of smooth curves 2d + 2g — 2, to 
degenerate curves. 

Also, some points of moduli spaces may have extra symmetries, and the moduli 
spaces are orbifolds (one needs to quotient by the symmetries), and intersection num- 
bers are counted with the symmetry factor, so that they can be rational numbers 
instead of integers. 

If we also fix the images of points Pi G C to be G P 1 , the Gromov-Witten 
invariants of P 1 are defined as the intersection numbers [59J: 



Those intersection numbers count a finite number of intersection points in 
Ji4 g>n (F l ,d), and thus they count certain curves C with marked points pi, ...,p n 
mapped to qx, . . . , q n . Those curves are instantons which extremize a topological sigma 




c 1 {C 1 ) k K..c 1 {C n f 



(4.131) 
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model action (type A topological string theory). To construct a curve which passes 
through given points, it is sufficient to know the first few terms in the Taylor expansion 
near those points, and thus the degree of tangency at those points, i.e. the number of 
vanishing coefficients in the Taylor expansion. Instantons are entirely encoded by those 
degrees, which form a partition of the total degree d, and the intersection numbers can 
be rewritten in terms of sums over partitions. 

In [50] H9], it was proved that the Gromov-Witten invariants of P 1 are: 

<T kl ...T km > d = YYJT-—T7j E V ^ H C *+lW ( 4 - 132 ) 
1 U\"* i" *■)■ x,\X\=d i 

This sum is nearly of the type of what we computed in this article, we only have to 
take derivatives with respect to the t k s at t k = 0. The general method for computing 
the derivatives of F g with respect to any parameter of the spectral curve is explained 
in [20], and can be applied here. It would be interesting to study in further details, the 
consequences of our formula for Gromov-Witten invariants of P 1 . 

4.4.2 Gromov-Witten invariants of CY 3-folds and topological strings 

See [43j [57J HQl [7J 08j [9] for an introduction to topological string theory and Toric 
Calabi-Yau 3-folds. 

For any p G Z, consider the toric Calabi Yau 3-fold X p , which was investigated by 

[21 HH SH H21 EZl EEJ: 

0(p-2)®0(-p) — >P X (4.133) 
which is a rank 2 bundle over P 1 . It is often represented by the toric fan diagram [57] : 




The partition function of the A-sigma model topological string theory in target space 
X p , is defined as the generating function of Gromov-Witten invariants of X p , i.e. a 
generating function for counting instantons, i.e. for computing intersection numbers: 

oo oo oo 

In {Z Xp {g s , t)) = 9 2 s 9 ~ 2 F g (t) = £ E *U*p) ^ ds 9 ' 2 (4-134) 

9=0 g=0 d=l 

Ng t d(X p ) are the Gromov Witten invariants of X p , i.e. in some sense the number of 
instantons (rational curves embedded in X p ) of genus g and degree d (the numbers 
may be rational because of orbifold points). 
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It can be computed with the so-called Topological Vertex method [3J, and it can be 
written as a sum of partitions with Plancherel measure [HI HU, [12] : 

Zx p (9s,t) = 52V 9 (\) ? HW)/, e^ , g = e"* (4.135) 

A 

where |A| = ^ Aj is the weight of the representation A, and C^A) = Aj(Aj — 2i + l) 
is the second Casimir, and V q (X) is the g-deformed Plancherel measure. 

Through mirror symmetry [501 El ED E3 EH], the type A topological string theory 
on X p , should be dual to a type B topological string theory on a mirror CY 3-fold, 
which we note X p , and whose equation is expected to be of the form: 

H(e*,e v )=Z{ (4.136) 

(with x,y,^X £ ^ 4 > so that this is locally a 3-dimensional complex submanifold of 
C 4 ), and where H(u,v) is a polynomial in both variables. Expressions of H(u,v) (up 
to symplectic transformations) were derived in [T2"| BTT]. 

F was computed in [121 ED EH], and Fi was computed in [UJ. Marcos Marino 
conjectured that higher F 5 's with g > 2 are given by the symplectic invariants of [20] 
for the curve H(e x ,e y ) = 0. 

Here, we have proved that conjecture for X p . 

To summarize, we have rewritten the type A model partition function of X p as 
a matrix integral, and therefore we obtained the topological small g s expansion from 
loop equations (in other words from Virasoro constraints). To leading order, the loop 
equations define a spectral curve (of mirror type). Once the spectral curve is known, 
all the topological expansion of the matrix integral, and thus the Gromov-Witten 
invariants, are then obtained from the method of [20J, as was conjectured by Marcos 
Marino [4T] . 

Finally, we check that the spectral curve of the matrix model (which coincides with 
the one found by [12] is symplecticaly equivalent to the mirror curve H(e x ,e y ) = 0, 
and therefore it indeed has the same F s 's. 

5 Conclusion 

In this article, we have shown that the topological expansion of a sum over partitions 
with the Plancherel measure: 

Z N (q;t k )= V(X) e~^ tJ ^ c ^ (5.137) 

n(X)<N 
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is a matrix integral, and has a topological expansion of the form: 



\nZ N {q-t k ) ~ J2^~ 9 F ^ ( 5 ' 138 ) 

9=0 



where the _F 9 's are the symplectic invariants of the spectral curve: 



x(z) = e~ u ° (z + -- ui) 



y(z) = \n(z) + \ J2i=iMz k - 



(5.139) 



Moreover, the same method also gives the large q asymptotic expansion, to all orders, 
of the density correlation functions. 

Our results were derived in a 1-cut regime, but it seems easy to extend them to 
any other regime, the only difference is that we would need ^-functions instead of the 
Zhukovski rational parametrization. 

We have also obtained similar results for the q— deformed Plancherel measure: 

Z N (q;t k )= V q (X)e-^^ tj ^ CkW , g = e"* (5.140) 

n(X)<N 

which is also a matrix model, whose topological expansion is of the form: 

oo 

\nZ N ~Y,g 2 s 9 ~ 2 F a (t k ) (5.141) 

9=0 

When ti = t,t 2 = p — 1, the sum is the generating function of Gromov-Witten 
invariants of the calabi-Yau 3-fold X p . The spectral curve is: 

x(z) = "-A>"-A> . . 

1 ' (1+ *» 2 , ^ , e- = 1 (1 - i)^ 2) (5-142) 



l/zz 



We have thus proved that can be obtained as a matrix integral (a unitary 
matrix integral, because the eigenvalues are constrained on a circle), and that the F g 's 
are the symplectic invariants of the mirror spectral curve H(e x ,e v ) = 0. This proves 
the conjecture of M. Marino [JT] for X p . 

Dijkgraaf and Vafa recentely noticed that that the symplectic invariants of [20J are 
also obtained from the Kodaira-Spencer theory [T7] , and thus we have proved that the 
type B topological string theory on the mirror of X p , is equivalent to a matrix model, 
and to the Kodaira-Spencer theory. 
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To go further in the same direction for all toric Calabi-Yau 3-folds, one would need 
to extend the method to plane partitions, because the topological vertex [3] is expressed 
as a sum of plane partitions, and it is only for X p that it reduces to usual partitions. 
Then one needs to "glue" topological vertices together. This is a work in progress [22], 
but with more involved tools because the spectral curve is no longer hyperellitical. 

Many other extensions of the present method can be explored. First, we have not 
studied in details the consequences of our method for crystal growth, longest increasing 
sequence, or T.A.S.E.P.. 

Also, we have not studied in detail the expectation values of Casimirs. That could 
have interesting conseqences in the Hurwitz numbers generating functions. Finaly, 
we have not computed the spectral curve of the q— defomed Plancherel measure, with 
higher Casimirs. The method works in the same way, but the solution is slightly more 
difficult to write. It would be interesting to continue this computation. 

The same kind of sums over partitions also appears in Seiberg-Witten theory, as 
discussed in |42j. We indeed find the same spectral curve and same prepotential F . 
The authors of [42] found their results from the integrable hierarchy structure, and 
indeed the symplectic invariants of |20j have such an integrable structure (see section 
9 in [20J). It would be interesting to explore further the computations for that case. 

Also, we have studied only 1-cut cases, which are the most relevant for combinatorics 
formal generating series. However, for many applications to statistical physics, we have 
convergent series, and it would be interesting to study phase transitions. When there 
are several cuts, the matrix integral no longer has a topological expansion, it receives 
an extra oscillatory part, as dicussed in [H 124"] . 
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